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ABSTRACT 

The rare isotope Li is made only by cosmic rays, dominantly in aa — > Li fusion 
reactions with ISM helium. Consequently, this nuclide provides a unique diagnostic of 
the history of cosmic rays in our Galaxy. The same hadronic cosmic-ray interactions 
also produce high-energy 7 rays (mostly via pp — > tt° — > 77). Thus, hadronic 7-rays 
and 6 Li are intimately linked. Specifically, 6 Li directly encodes the local cosmic-ray 
' fluence over cosmic time, while extragalactic hadronic 7 rays encode an average cosmic- 

ray fluence over lines of sight out to the horizon. We examine this link and show 
how 6 Li and 7-rays can be used together to place important model-independent limits 
on the cosmic-ray history of our Galaxy and the universe. We first constrain 7-ray 
production from ordinary Galactic cosmic rays, using the local 6 Li abundance. We 



> 



Or 



find that the solar 6 Li abundance demands an accompanying extragalactic pionic 7-ray 
intensity which exceeds that of the entire observed EGRB by a factor of 2 — 6. Possible 
explanations for this discrepancy are discussed. We then constrain Li production using 
recent determinations of extragalactic 7-ray background (EGRB). We note that cosmic 
^ ■ rays created during cosmic structure formation would lead to pre-Galactic Li production, 

which would act as a "contaminant" to the primordial 7 Li content of metal-poor halo 
stars; the EGRB can place an upper limit on this contamination if we attribute the 
entire EGRB pionic contribution to structure forming cosmic rays. Unfortunately, the 
uncertainties in the determination of the EGRB are so large that the present 7-ray data 
cannot guarantee that the pre-Galactic Li is small compared to primordial 7 Li; thus, 
an improved determination of the EGRB will shed important new light on this issue. 
Our limits and their more mo del- dependent extensions will improve significantly with 
additional observations of 6 Li in halo stars, and with improved measurements of the 
EGRB spectrum by GLAST. 

Subject headings: cosmic rays - gamma rays - nuclear reactions, nucleosynthesis, abun- 
dances 
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1. Introduction 

The origin and history of cosmic rays has been a subject of intensifying interest. For more 
than a decade, a large body of work has focused on the light elements Li, Be, and B (LiBeB) as 
signatures of cosmic-ray interactions with the diffuse gas (for a recent review see Casse, Vangioni- 
Flam, & Audouze 2001). LiBeB abundances in Galactic halo stars have been used to probe the 
history of cosmic rays in the (proto-) Galaxy More recently, a great deal of attention has been 
focused on high-energy 7-rays also produced in interactions during cosmic-ray propagation. Here, 
we draw attention to the tight connection between these observables, particularly between 7-rays 
and 6 Li. 

The abundances of LiBeB nuclei encode the history of cosmic ray exposure in local matter. In 
the past 15 years or so, measurements of LiBeB in the Sun and in Galactic disk have been joined 
by LiBeB observations in halo stars; these offer particularly valuable information about cosmic-ray 
origins and interactions in Galactic and proto-Galactic matter. In particular, different scenarios 
for cosmic ray origin lead to different LiBeB trends, which have been modeled and compared 
with observations (see, e.g., Vangioni-Flam & Casse 2001; Fields & Olive 1999a; Ramaty, Scully, 
Lingenfelter, & Kozlovsky 2000, and references therein). For the purposes of this paper, the details 
of these models are less important than the following basic distinction: all LiBeB species are 
produced as cosmic rays interact with interstellar gas and fragment- "spall" -heavy nuclei, e.g., 
p + O — > 9 Be. However, the fusion processes a + a — > 6 ' 7 Li yield lithium isotopes exclusively, and 
indeed dominate the cosmic-ray production of Li (Steigman & Walker 1992; Montmerle 1977c). This 
makes cosmic-ray lithium production particularly "clean" since its evolution depends uniquely on 
its exposure to cosmic rays, and unlike Be and B, does not depend on the ambient heavy element 
abundances. 

Cosmic-ray interactions provide the only known source for the nucleosynthesis of 6 Li, 9 Be, and 
10 B, making these species ideal observables of cosmic ray activity. 1 The story is more complex 
for n B, which can also be produced in core-collapse supernovae by the "neutrino process" (e.g., 
Woosley, Hartmann, Hoffman, & Haxton 1990; Yoshida, Terasawa, Kajino, & Sumiyoshi 2004). 
Finally, 7 Li has the most diverse lineage. In the early Galaxy, and hence in halo stars, 7 Li is 
dominated by the contribution from primordial nucleosynthesis (e.g., Cyburt, Fields, & Olive 2003, 
and references therein), with a small contribution from cosmic-ray fusion as well as the neutrino 
process (Ryan et al. 2000). At late times, and hence in disk stars including the Sun, 7 Li has 
important and probably dominant contributions from longer-lived, low-mass stars (though the 
specific site remains controversial: Romano, Matteucci, Ventura, & D'Antona 2001; Travaglio et 
al. 2001). In this paper we will build on the work of Suzuki & Inoue (2002) to point out the 



In fact, a pre-Galactic component of Li can be produced in some scenarios in which dark matter decays via 
hadronic (Dimopoulos, Esmailzadeh, Starkman, & Hall 1988) or electromagnetic (Jedamzik 2000; Kawasaki, Kohri, 
& Moroi 2001; Cyburt, Ellis, Fields, & Olive 2003) channels. Such scenarios are constrained via their effects on the 
other light elements, but some level of 6 Li production is hard to rule out completely. 
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possible importance of another, pre- Galactic, source of cosmic-ray 7 Li and 6 Li, which could confound 
attempts to identify the pre-Galactic Li abundance with the primordial component. We cannot 
rule out (or in!) this possible source, but we will constrain it using observations of 7-rays. 

Cosmic-ray interactions with interstellar gas produce not only LiBeB, but also inevitably pro- 
duce 7-rays. Cosmic rays in the Galactic disk today lead to pronounced emission seen in the 
Galactic plane (Hunter et al. 1997). Cosmic ray populations in (and between!) external galaxies 
would contribute to a diffuse extragalactic 7-ray background (hereafter the EGRB). The existence 
of an EGRB was already claimed by some of the first 7-ray observations (Fichtel, Kniffen, & Hart- 
man 1973). The most recent high-energy (i.e., roughly in the 30 MeV - 30 GeV range) 7-ray 
observations are those of the EGRET experiment on the Compton Gamma-Ray Observatory, and 
the EGRET team also found evidence for a EGRB (Sreekumar et al. 1998). The intensity, energy 
spectrum, and even the existence of an EGRB are not trivial to measure, as this information only 
arises as the residual after subtracting the dominant Galactic foreground from the observed 7-ray 
sky. The procedure for foreground subtraction is thus crucial, and different procedures starting 
with the same EGRET data have arrived at an EGRB with a lower intensity and different spec- 
trum (Strong, Moskalenko, & Reimer 2004), or have even failed to find evidence for an EGRB at 
all (Keshet, Waxman, & Loeb 2003). Despite these uncertainties, we will see that the EGRB (or 
limits to it) and Li abundances are mutually very constraining. 

Whether or not an EGRB has yet been detected, at some level it certainly should exist. 
EGRET detections of individual active galactic nuclei (blazars) as well as the Milky Way and the 
LMC together guarantee that unresolved blazars (e.g., Stecker & Salamon 1996; Mukherjee & Chi- 
ang 1999), and to a lesser extent normal galaxies (Pavlidou & Fields 2002), will generate a signal at 
or near the levels claimed for the EGRB. Many other EGRB sources have been proposed, but one 
of the promising has been a subject of intense interest recently: namely, 7-rays originating from a 
cosmological component of cosmic rays. This as- yet putative cosmic-ray population would originate 
in shocks (Miniati et al. 2000; Keshet et al. 2003; Ryu, Kang, Hallman, & Jones 2003) associated 
with baryonic infall and merger events during the growth of large-scale cosmic structures. Diffusive 
shock acceleration (e.g., Blasi 2004; Kang, Jones, & Gieseler 2002; Jones k, Ellison 1991) would then 
generate a population of relativistic ions and electrons. Gamma-ray emission would then follow 
from inverse Compton scattering of electrons off of the ambient photon backgrounds and from tt° 
production due to hadronic collisions (Loeb & Waxman 2000). The most recent semi-analytical 
and numerical calculations (Gabici & Blasi 2003; Miniati 2002) suggest that this "structure form- 
ing" component to the EGRB is likely below the blazar contribution, but the observational and 
theoretical uncertainties here remain large; upcoming 7-ray observations by GLAST (Gehrels & 
Michelson 1999) will shed welcome new light on this problem. 

The link between the nucleosynthesis and 7-ray signatures of cosmic-ray history has been 
pointed out by others in multiple contexts. We note in particular the prescient work of Montmerle 
(1977a, b,c), who in a series of papers considered the implications of a hypothetical population of 
"cosmological cosmic rays" in addition to the usual Galactic cosmic rays. Montmerle's analysis 
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is impressive in its foresight and its breadth. Montmerle (1977a) develops the formalism for a 
homogeneous population of cosmological cosmic rays (assumed to be created instantaneously at 
some redshift), and describes their propagation in an expanding universe, as well as their light 
element and 7-ray production. He identifies the tight connection between 6 Li and extragalactic 
7-rays, and exploits this connection to use the available EGRB data to constrain Li production 
for a variety of different assumptions. A particularly pertinent case involves an EGRB near the 
levels discussed today ("normalization 2" in Montmerle's parlance), coupled with a cosmic baryon 
density close to modern values (e.g., Spergel et al. 2003; Cyburt, Fields, & Olive 2003). Under these 
conditions, Montmerle (1977b) finds that cosmological cosmic ray activity sufficient to explain the 
EGRB leads to a present 6 Li abundance that is about an order of magnitude smaller than the 
solar abundance. This result foreshadows an important conclusion we will find: if the solar 6 Li 
abundance is produced by Galactic cosmic rays, then the associated pionic 7-ray production exceeds 
the entire EGRB by at least a factor of 2. 

More recently, studies of structure formation cosmic rays have focused primarily on their 7- 
ray signatures. However, recently Suzuki & Inoue (2002) also proposed using 6 Li as a diagnostic 
of shock activity in the Local Group. These authors note that the resulting 6 Li abundances in 
halo stars could be used to probe the shocks and resulting cosmic rays in proto-Galactic matter. 
We also will draw on this idea, with an emphasis on the fact that pre-Galactic Li production 
would be (by itself) observationally indistinguishable from the primordial 7 Li production from big 
bang nucleosynthesis. Thus we attempt to use EGRB data to constrain this possibility, but find 
that current data is unable to rule out a significant contribution to halo star Li abundances from 
structure formation cosmic rays. 

Our work thus follows these pioneering efforts, further emphasizing and formally exploring the 
intimate connection between cosmic-ray nucleosynthesis and high-energy 7-ray astrophysics. In §2, 
we formally show and discuss the generality and tightness of the 6 Li-7 connection, and in §3 the 
relevant observations are reviewed. In §4 we use the theory of LiBeB production by Galactic cosmic 
rays to deduce the minimal contribution the EGRB. In §5 we exploit the link to use the observed 
EGRB to limit the SFCR contribution to pre-Galactic lithium. Discussion and conclusions appear 
in §6. 

2. The Gamma-Ray — Lithium Connection: Formalism 

Before doing a detailed calculation let us first establish a simple, back of the envelope, connec- 
tion between gamma-rays and lithium. We know that low energy (~ 10 — 100 MeV/nucleon) 
hadronic cosmic rays produce lithium through aa — > 6 ' 7 Li + • • •. But higher-energy (> 280 
MeV/nucleon) cosmic rays also produce 7-rays via neutral pion decay: pp — > 7r° — > 77. Because 
they share a common origin in hadronic cosmic ray interactions, we can directly relate cosmic ray 
lithium production to "pionic" gamma-rays. The cosmic-ray production rate of 6 Li per unit vol- 
ume is c/( 6 Li) = a aa ^6 Li & a n a , where $ a is the net cosmic ray He flux, n a is the interstellar He 
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abundance, and a^^e,^ is the cross section for 6 Li production, appropriately averaged over the 
cosmic-ray energy spectrum (detailed definitions and normalization conventions appear in Appendix 
A). Thus, the 6 Li mole fraction Y 6 = n 6 /n h is just Y( 6 Li) ~ f ||<76 Li ~ ya,cvY a ,is^ aa ^u^pto- 
where y a ^ = $ a /® p ps (He/H) ism . 

On the other hand, the cosmic-ray production rate of pionic 7-rays is just the pion production 
rate times a factor of 2, that is, = 2o pp ^ 7r o 3> P: 

,cr^-p,ism- Integrated over a line of sight towards the 
cosmic particle horizon, this gives a EGRB intensity l 1 ~ c f dtq^/A-K ~ 2cr pp ^ 7r oc$ p to- Thus we see 
that both the 6 Li abundance and the 7-ray intensity have a common factor of the (time-integrated) 
cosmic-ray flux, and so we can eliminate this factor and express each observable in terms of the 
other: 

27]- a 6Lj 

^6t j ~ ya,cr^a,gas— — pp ^7 (1) 

From eq. (1) we see that the connection between cosmic-ray lithium production and pionic gamma- 
ray flux is straight forward. 

This rough argument shows the intimacy of the connection between 6 Li and pionic 7-rays. 
However, this simplistic treatment does not account for the expansion of the universe, nor for time- 
variations in the cosmic-cosmic ray flux, nor for the inhomogeneous distribution of sources within 
the universe. We now include these effects in a more rigorous treatment. 

For Li production at location x, the production rate per unit (physical) volume is 

qu(x) = (T aa <S>a{x)n a>gas (x) = y a , C rY^ sm a aa ^p(x)n h:gas (x) = /i(x)r Li (x)n b (x) (2) 

Here, y QiCr = (a/p) CI is the cosmic-ray He/H ratio, and is assumed to be constant in space and 
time. 2 The target density of (interstellar or intergalactic) helium is n aigas , which we write in terms 
of its ratio Y^ sm = nne/^b to the baryon density. We take Y^ sm ~ 0.06 to be constant in space and 
time, but we do not assume this for the baryon density nt,(x). The baryonic gas fraction 

At = ™ b ,gas/nb ( 3 ) 

accounts for the fact that not all baryons need to be in a diffuse form. Finally, we will find it 
convenient to write qu(x) in terms of the local baryon density and the local Li production rate 
Tu(x) per baryon. 



With these expressions, we have 



jYUx) = Kx)TUx) (4) 



2 That is, we ignore the small non-primordial 4 He production by stars, and we neglect any effects of H and He 
segregation. Both of these should be quite reasonable approximations. 
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which we can solve to get 

Y Li (x,t) = [ dt' ti(x,t') r Li (x,0 (5) 
J o 

= y a ,crYra aa [ dt' fi(x,t')^(x,t') (6) 
Jo 

= y*,crYi sm a aa F p (x,t) (7) 

where F p (x,t) = f^dt' /J-(x,t') 3>p*(x, t') is the local proton fluence (time-integrated flux), weighted 
by the gas fraction. Thus we see that Li (and particularly 6 Li) serves as a "cosmic-ray dosimeter" 
which measures the net local cosmic ray exposure. 

We now turn to 7 rays from hadronic sources, most of which come from neutral pion production 
and decay: pp — > 7r° — > 77. The extragalactic background due to these process is expected 
to be isotropic (at least to a good approximation). In this case, the total 7-ray intensity 7 7 = 
dNj/dAdtdQ, integrated over all energies, is given by an integral 

m = ^[dt> q c ™\t>) (8) 

of the sources over the line of sight to the horizon. We are interested in particular in the case of 
hadronic sources, so that q com = a 3 q is the total (energy-integrated) comoving rate of hadronic 
7-ray production per unit volume; here a is the usual cosmic scale factor, which we normalize to 
a present value of ao = a(to) = 1. A formal derivation of eq. (8) appears in Appendix B, but 
one can arrive at this result from elementary considerations. Namely, note that the comoving 
number density of photons produced at any point is just n 7)Com = J*g <? 7 ° m dt' . We neglect photon 
absorption and scattering processes, and thus particle number conservation along with homogeneity 
and isotropy together demand that the comoving number density of ambient photons at any point 
is the same as the comoving number density of photons produced there. Furthermore, the total 
(energy-integrated) photon intensity is also isotropic and thus by definition is I 1 = n 7jCom c/47r, 
which is precisely what we find in eq. (8). 

The comoving rate of pionic 7-ray production per unit volume at point s is 

q^(s, t) = a^ p (s, t)n%™ s (s, t) = fi(s, t)a^ p (s, t)n% m (s, t) (9) 

where nu is the (comoving) hydrogen density, and & p = 47r j I p (e Tes t)de res t is the total (integrated 
over rest-frame energy e rcst ) omnidirectional cosmic ray proton flux. The flux-averaged pionic 7-ray 
production cross section is 

_ oc > ot I de TCSt I p (e res t) CnO' 7T o(€ res t) . . 

07 = I^Qa^o = 2£ a -TT- — r (W) 

J "Crest JpV^rest ) 

where the factor of 2 counts the number of photons per pion decay, a w o is the cross section for 
pion production and C?r is the pion multiplicity, and the factor £ a = 1.45 accounts for pa and aa 
reactions (Dermer 1986). 
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Then we have 

fi com / — \ 

/ ( t ) = ^Y h <t 7 / A ,x(|) t) = ^a 7 yif m F p (t) (11) 

4vr Jo n bi0 47r 

where 

r-t com/^\ 

F p (t) = / dt fi(s) -«2-Ai d> p(s -) (12) 
Jo "b,0 

is a mean value of the cosmic-ray fluence along the line of sight, where the average is weighted by 
the gas fraction and the ratio n^ >m (s)/n^ Q of the local baryon density along the photon path. Note 
that the 7-ray sources are sensitive to the overlap of the cosmic-ray flux with the diffuse hydrogen 
gas density, and thus need not be homogeneous. Even so, we still assume the ERGB intensity to 
be isotropic, which corresponds to the assumption that the line-of-sight integral over the sources 
averages out any fluctuations. 

One further technical note: I 1 = i" 7 (> 0) = J °° de 7 i" 7 (e 7 ) represents the total pionic 7-ray flux, 
integrated over photon energies. While this quantity is well-defined theoretically, real observations 
have some energy cutoff, and thus report I 7 (> eo) = J e °° c/e 7 J 7 (e 7 ), typically with eo = 100 MeV. 
But the spectrum of pionic 7-rays will be shifted towards lower energies if they originate from a 
nonzero redshift. Thus it is clear that 7-ray intensity I 7 , integrated above some energy eo 7^ 0, will 
be redshift-dependent. A way to eliminate this z-dependence is to include all pionic 7-rays, that is 
to take 7 7 = 7 7 (> GeV), i.e., to take eo = 0. As discussed in more detail in Appendix B, the 6 Li-7 
proportionality is only exact for 7 7 (> 0), as this quantity removes photon redshifting effects which 
spoil the proportionality for eo 7^ 0. Thus we will have to use information on the pionic spectrum 
to translate between / 7 (> eo) and i" 7 (> 0); these issues are discussed further in §3.1. 

Thus we see that the lithium abundance and the pionic 7-ray intensity (spectrum integrated 
from energy) arise from very similar integrals, which we can express via the ratio 

J 7 (t) _ n b c 07 F p (t) ^ 



Yi{x,t) 47ry Q , !Cr y ai ; sm cr % aoi F p (x,t) 

where i denotes 6 Li or 7 Li. Note that this "7-to-lithium" ratio has its only significant space and 
time dependence via the ratio F p (t)/F p (x, t) of the line-of-sight baryon-averaged fluence to the local 
fluence. 3 

The relationship expressed in eq. (13) is the main result of this paper, and we will bring this 
tool to bear on Li and 7-ray observations, using each to constrain the other. To do this, it will be 
convenient to write eq. (13) in the form 

Yj(x,t) F p (t) 

w- 1 * v hQ Fp m (14) 



3 In fact, the ratio also depends on the shape of the cosmic-ray spectrum (assumed universal), which determines 
the ratio of cross sections. We will take this into account below when we consider different cosmic-ray populations. 
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where the scaling factor 

h,i = j — Y it Q (15) 

is independent of time and space, and only depends, via the ratio of cross sections, on the shape of 
the cosmic-ray population considered. Table 1 presents the values of io,i for the different spectra 
that will be considered in the following sections. Values of the scaling factor were obtained by 
using photon multiplicity £ 7 = 2, £ Q = 1.45, baryon number density = 2.52 x 1CP 7 cm" 3 , CR 
and ISM helium abundances = = 0.1 and solar abundances = 1-53 x 10~ 10 and 

y7"Li = 1-89 x 10~ 9 (Anders & Grevesse 1989). For the 7r° and lithium production cross-sections, 
we used the fits taken from Dermer (1986) and Mercer et al. (2001), and from that obtained the 
ratios of flux-averaged cross-sections for different spectra, and these are also presented in Table 1. 



Table 1: Lithium and 7-ray Scalings and Production Ratios 



Cosmic-Ray 


Io,6 


I0.7 








Population 


[cm~ 2 s~ 




aa IVV 


TV 


7 Li/ 6 Li 


GCR 


9.06 x 10~ 5 


8.36 x 10~ 4 


0.21 


0.28 


1.3 


SFCR 


1.86 x 10~ 5 


1.15 x 10~ 4 


1.02 


2.03 


2.0 



Table 1 shows that the different cosmic-ray spectra lead to very different Li-to-7 ratios. For 
example, the 6 Li-to-7 ratio o-Q^./a pp is almost a factor of 5 higher in the SFCR case than in the 
GCR case. The reason for this stems from the different threshold behaviors and energy dependences 
of the Li and ir° production cross sections. Li production via aa fusion has a threshold around 10 
MeV/nucleon, above which the cross section rapidly rises through some resonant peaks. Then be- 
yond ~ 15 MeV/nucleon, the cross section for 6 Li drops exponentially as e~ E / 16 MeV/nucieon^-y ercer 
et al. 2001), rapidly suppressing the importance of any projectiles with E 3> 16 MeV/nucleon. 
Thus, as has been widely discussed, Li production is a low-energy phenomenon for which the 
important projectile energy range is roughly 10 — 70 MeV/nucleon. 

On the other hand, pp — > it production has a higher threshold of 280 MeV, and the effective 
cross section C n o-p p rises with energy up to and beyond 1 GeV. Neutral pion production is thus a 
significantly higher-energy phenomenon. 

These different cross section behaviors are sensitive to the differences in the two cosmic-ray 
spectra we adopt. On the one hand, we adopt a GCR spectrum that is a power law in total energy: 
4> P (E) oc (m p + E)~ 2 ' 75 , a good approximation to the locally observed (i.e., propagated) spectrum. 
This spectrum is roughly constant for E < m p . Thus, there is no reduction in cosmic-ray flux 
between the Li and 7r° thresholds. Furthermore the flux only begins to drop far above the ir° 
threshold at 280 MeV, so that there is significant pion production over a large range of energies, 
in contrast to the intrinsically narrow energy window for Li production. As a result of the effects, 
o-^./a^ < 1 for the GCR case. 
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In contrast, the SFCR flux is taken to be the standard result for diffusive acceleration due 
to a strong shock: namely, a power law in momentum <p(E) oc p(E)~ 2 . This goes to <f> oc E^ 1 at 
E ^ m p , and </> oc E~ 2 at higher energies. This spectrum thus drops by a factor of 28 between the 
Li and tt° thresholds, and continues to drop above the 7r° threshold, offsetting the rise in the pion 
cross section. This behavior thus suppresses 7r° production relative to the GCR case, and thus we 
have a significantly higher <Tg£./cr?? p ratio. As we will see, these ratios-and the differences between 
them-will be critical in deriving quantitative constraints. 

3. Observational Inputs 

We have seen that the EGRB intensity and lithium abundances are closely linked. Here we 
collect information on both observables. 



3.1. The Observed Gamma-Ray Background and Limits to the Pionic Contribution 

Ever since 7-rays were first observed towards the Galactic poles as well as in the plane (Fichtel, 
Kniffen, &: Hartman 1973), the existence of emission at high Galactic latitudes has been regarded as 
an indication of an EGRB. However, any information regarding the intensity, energy spectrum, and 
even the existence of the EGRB is only as reliable as the procedure for subtracting the Galactic 
foreground. Such procedures are unfortunately non-trivial and model-dependent. The EGRET 
team (Sreekumar et al. 1998) used an empirical model for tracers of Galactic hydrogen and starlight, 
and found evidence for an EGRB which dominates polar emission. Other groups have recently 
presented new analyses of the EGRET data. In a semi-empirical approach using a model of Galactic 
7-ray sources, Strong, Moskalenko, & Reimer (2004) also find evidence for an EGRB, but with a 
different energy spectrum and a generally lower intensity than the Sreekumar et al. (1998) result. 
Finally, Keshet, Waxman, & Loeb (2003) find that the Galactic foreground is sufficiently uncertain 
that its contribution to the polar emission can be significant, possibly saturating the observations. 
Consequently, the Keshet, Waxman, & Loeb (2003) analysis is unable to confirm the existence of 
an EGRB in the EGRET data; instead, they can only to place upper limits on the EGRB intensity. 

It was recently shown by Prodanovic & Fields (2004) that a model-independent limit on the 
fraction of EGRB flux that is of pionic origin (gamma rays that originate from ir° decay) can be 
placed. Their limit comes from noticing that the EGRB shows no strong evidence of the distinctive 
shape pionic 7-ray spectral peak at m n o/2, the "pion bump." Thus by comparing the shapes of the 
observed EGRB and theoretical pionic gamma-ray spectrum, they were able to maximize the pionic 
flux so that it stays below the observed one. This procedure allowed them to place constraints on 
the maximal fraction of EGRB that can be of pionic origin. 

For the pionic 7-ray source-function, Prodanovic & Fields (2004) used a semi-analytic fit from 
the Pfrommer & Enfilin (2003) paper and the Dermer (1986) model for the production cross section. 
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A key feature of the pionic 7-ray spectrum is that it approaches a power law at both high and low 
energies, going to e a ~< for e <C and to e - " 7 for e 3> m n /2. In Dermer's model, the 7-ray 

spectral index a 7 is equal to the cosmic-ray spectral index. Prodanovic & Fields (2004) adopted 
the value a 7 = 2.2 for pionic extragalactic 7-rays, which is consistent with blazars and structure- 
forming cosmic rays as their origin. In this simple analysis, Prodanovic & Fields (2004) used a 
single-redshift approximation, that is, they assumed that these 7-rays are all coming from one 
redshift, and thus their limit on the maximal pionic fraction is a function of z. 

To obtain the EGRB spectrum from EGRET data, a careful subtraction of Galactic foreground 
is needed. Prodanovic & Fields (2004) considered two different EGRB spectra and obtained the 
following limit : for the Sreekumar et al. (1998) spectrum they found that the pionic fraction 
of the EGRB (integrated spectra above 100 MeV) can be as low as about 40% for cosmic rays 
that originated at present, to about 90 % for z = 10; for the more shallow spectrum of Strong, 
Moskalenko, h Reimer (2004) they found that pionic fraction can go from about 30% for z = up 
to about 70% for z = 10. However, the Keshet, Waxman, & Loeb (2003) analysis of the EGRET 
data implies that the Galactic foreground dominates the 7-ray sky so that only an upper limit on 
the EGRB can be placed, namely i" 7 (> lOOMeV) < 0.5 x 10 -5 cm~ 2 s _1 sr _1 . Thus, in this case, 
we were not able to obtain the pionic fraction. 

However, to be able to connect the pionic 7-ray intensity I 7 with lithium mole fraction Y{ as 
shown in (13), J 7 must include all of the pionic 7-rays, that is, the spectrum has to be integrated 
from energy eo = 0. The upper limit to the pionic 7-ray intensity above energy eo for a given 
redshift can be written as 



where f n (> eo,z) is the upper limit to the fraction of pionic 7-rays (Prodanovic & Fields 2004), 
I 7 bs (> eo) is the observed intensity above some energy, while tp[e(l + z)\ is the semi-analytic fit 
for pionic 7-ray spectrum (Pfrommer & Enfilin 2003) which is maximized with A/" max normalization 
constant. An upper limit to the pionic 7-ray intensity that covers all energies J 7 (> 0, z), follows 
immediately from the above equations: 



7 7 (> e ) 



)C S (> *>) 



(16) 




(17) 



i" 7 (> 0,z) = f w (> e ,z 



)C S (> eo) 



f Q y>[e(l + z)]de 
Lp[e(l + z)]de 



(18) 



Now this is something that is semi-observational and can be easily obtained from 7-ray intensity 
observed above some energy, and from Prodanovic & Fields (2004) and Pfrommer & Enfilin (2003) 
results. 
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3.2. The Primordial Lithium Abundance 

Given the EGRB intensity, we will infer the amount of associated lithium production. It will be 
of interest to compare this not only to the solar abundance, but also to the primordial abundance 
of 7 Li. Metal-poor halo stars (extreme Population II) serve as a "fossil record" of pre-Galactic 
lithium. Ryan et al. (2000) find a pre-Galactic abundance 

^) = (1-23^) x 10- 10 (19) 

/ pre— Gal,obs 

based on an analysis of very metal-poor halo stars. On the other hand, one can use the WMAP 
(Spergel et al. 2003) baryon density and BBN to predict a "theoretical" (or "CMB-based") primor- 
dial 7 Li abundance (Cyburt, Fields, & Olive 2003): 

C-^) = (3.82±g;g) x 10- 10 (20) 

V n J BBN.thy 

These abundances are clearly inconsistent. Possible explanations for this discrepancy include 
unknown or underestimated systematic errors in theory and/or observations or new physics; these 
are discussed thoroughly elsewhere (see, e.g., Cyburt, Fields, & Olive 2003, and refs therein). For 
our purposes, we will acknowledge this discrepancy by comparing pre-Galactic lithium production 
by cosmic rays with both the observed and CMB-based Li abundances. 



4. 6 Li and Gamma-Rays From Galactic Cosmic Rays 

We have shown that 6 Li abundances and extragalactic 7-rays are linked because both sample 
cosmic-ray fluence, and now apply this formalism to 7-ray and 6 Li data. In this section we turn 
to the hadronic products of Galactic cosmic rays, which are believed to be the dominant source of 
6 Li, but a sub-dominant contribution to the EGRB. 



4.1. Solar 6 Li and Gamma-rays 

We place upper limits on the lithium component of GCR origin by using the formalism es- 
tablished in earlier sections. To be able to find I^/Ye^ from eq. (13) we assume that ratio of 
cosmic-ray fluence along the line of sight (weighted by gas fraction) to the local cosmic-ray fluence 
is F p (t)/F p (x,t) pa 1. That is, we assume that the Milky Way fluence is typical of star forming 
galaxies, i.e., that the 7-luminosities are comparable: Lmw ~ (L) ga \- Note that in the most simple 
case of a uniform approximation (cosmic-ray flux and gas fraction the same in all galaxies), the 
two fluences would indeed be exactly equal. 
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Taking the solar 6 Li abundance and ■) / (&w P ) = 0.21 for the ratio of GCR flux averaged 



'l.i 

cross-sections, we can now use eq. (14) to say that / 7 7r o(e > 0) = 9.06 x 10 _5 cm~ 2 s _1 sr _1 is the 
hadronic 7-ray intensity that is required if all of the solar 6 Li is made via Galactic cosmic-rays. 

We wish to compare this 6 Li-based pionic 7-ray flux to the observed EGRB intensity I° bs (e > 
eo). However, eq. (14) gives the hadronic 7-ray intensity integrated over all energies, whereas the 
observed one is above some finite energy. Thus we have to compute 

f del e n 

/ 7i7r o(e>e ) = J 7|7r o(e>0) 7 ' (21) 
= 9.06 x l(T 5 cm- 2 s" 1 sr' 1 ^" f^'* (22) 

Jo ""e,7r 

We follow the model of Pavlidou & Fields (2002) to calculate the GCR emissivity over the history 
of the universe. The source function q^ om (equivalent to eq. 9) is given by a coarse-graining over 
galactic scales, so that 

q^z,e)=L,(e)n™?(z) (23) 

where L 7 is the average galactic 7-ray luminosity (by photon number), and n c °™(z) is the mean 
comoving number density of galaxies. The key assumptions for the luminosity L 7 are: (1) that 
supernova explosions provide the engines powering cosmic-ray acceleration, so that the cosmic-ray 
flux $ oc Y> scales with the supernova rate and thus the star formation rate tp', (2) that the targets 
come from the gas mass which evolves following the "closed box" prescription; and (3) that the 
Milky Way luminosity represents that of an average galaxy. With these assumptions we have that 
L 7 oc fitp, and thus that q^ om oc fj,p*, where is the cosmic star formation rate. 

Following Pavlidou & Fields (2002), the specific form of I en is expressed in terms of the 
present day Milky Way gas mass fraction //o,mw> cosmic star- formation rate p*(z), Milky Way 
gamma-ray (number) luminosity L 7i mw(z, E), cosmology Q\ and Q m , and integrated up to z*, the 
assumed starting redshift for star formation. For this calculation we adopt the following values: 
Mo,mw = 0.14, f^A = 0.7, Q m = 0.3 and = 5. For the cosmic star formation rate we use the 
dust-corrected analytic fit from Cole et al. (2001). Finally we need the (number) luminosity of 
pionic gamma-rays which we can write as 

L 7i mw(^, E) = T^Np = ^N p cx $M gas (24) 

Tip 

where n p is the proton number density in the Galaxy, N p is the total number of protons in the 
Galaxy, while g 7i7r [s _1 GeV~ 1 cm~ 3 sr~ 1 ] is the source function of gamma-rays that originate from 
pion decay adopted from Pfrommer & Enfilin (2003). Notice that in equation (21) we have the 
ratio of two integrals where integrands are identical, thus normalizations and constants will cancel 
out. Therefore, instead of using the complete form of L 7) mw(2,^) we need only use the spectral 
shape of the pionic gamma-ray source function (Pfrommer & Enfilin 2003), that is, only the part 
that is energy-and redshift-dependent. 
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Finally then, we find 

7 7i7r o(e > O.lGeV) = 3.22 x 10~ 5 cm" 2 s" 1 sr" 1 (25) 

which we can now compare to the observed EGRB values I° bs (e > O.lGeV) that are given in the 
first column of Table 2. As one can see, our pionic EGRB gamma-ray intensity is between 2 and 6 
times larger than the entire observed value! 

We thus conclude that the solar 6 Li abundance, if made by GCRs as usually assumed, seems 
to demand an enormous diffuse pionic 7-ray contribution, far above the entire EGRB level. How 
might this discrepancy be resolved? One explanation follows by dropping our assumption that 
-F(xmWj^o) = ^avg(^o)) i- e -> that the baryon-weighted Milky Way GCR fluence is the same as the 
cosmic mean for star-forming galaxies. Note that we have F = J dt/i<& oc J dt (ipM gas ), where tp is 
the global galactic star formation rate (assuming $ cx ip), and M gas the galactic gas content. If our 
Galaxy has an above-average star formation rate and/or gas mass, this will increase the local 6 Li 
production relative to the average over all galactic populations, and thus lead to an overestimate 
of the EGRB. 

In this connection it is noteworthy to compare our 6 Li-based estimate of the galactic EGRB 
contribution to the work of Pavlidou & Fields (2002). That calculation adopted the same model 
for the redshift history of cosmic-ray flux and interstellar gas, and so only differed from the present 
calculation in the normalization to Galactic values. Pavlidou & Fields (2002) normalized to the 
present Galactic 7-ray luminosity. This amounts to a calibration not to the time-integrated cosmic- 
ray fluence, but rather to the instantaneous cosmic ray flux, as determined by the Dermer (1986) 
emissivity, a Galactic gas mass of 1O 1O M , and an estimate of the present Galactic star formation 
rate. This normalization gave a galactic EGRB component which at all energies lies below the total 
(Sreekumar et al. 1998) background. Our calculation is normalized to solar 6 Li, which is a direct 
measure of Galactic (or at least solar neighborhood) cosmic-ray fluence, and which contains fewer 
uncertainties than the factors entering in the Pavlidou & Fields (2002) result. Yet surprisingly, the 
6 Li-based fluence result gives a high pionic EGRB, while the more uncertain normalization gives 
an acceptable result. 

Can we independently test whether our Galaxy has had an above-average cosmic-ray exposure? 
This present a challenge, as we require an integral measure of cosmic-ray activity, which is readily 
available locally but difficult to obtain in external galaxies. The best candidates are the LiBcB 
isotopes; 6 Li is ideal for the reasons we have outlined, but is not accessible in stars bright enough 
to be see in even the nearest external galaxies. The best hope then would be for measurements of 
9 Be, in the Local Group or beyond. In the SMC, such measurements have placed interesting limits 
on boron abundances Brooks et al. (2002), though the presence of the neutrino-process production 
of n B makes boron observations more difficult to interpret than beryllium. 

Another explanation for the high intensity of eq. (25) stems from noting that the required 
ERGB intensity is much larger for GCRs than that one would infer from SFCRs, due to the large 
difference in the 7r°/ 6 Li ratio for these two spectra. Were the Milky Way spectra is atypically skewed 
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to high energies, we would overestimate the 7r°/ 6 Li ratio and thus the EGRB contribution. This 
possibility (which we regard as less likely than the previous one) could be tested by observations of 
cosmic-ray spectra in external Galaxies, e.g., by 7-ray observations of Local Group galaxies such 
as GLAST should perform Pavlidou & Fields (2001). 

Finally, a related but more unconventional view would be that 6 Li is in fact primarily made 
by SFCRs themselves, rather than by GCRs. This suggestion is further discussed and constrained 
below, §5. 

We close this subsection by noting that if the depletion of 6 Li is taken into consideration, one 
might use 6 Li abundance larger than solar. In that case one would find that the accompanying 
pionic EGRB gamma-ray intensity is more than 2-6 times greater than the observed EGRB. 



4.2. The Observed EGRB and Non-Primordial 6 Li 

We can exploit equation (14) in both directions. Here we use the observed EGRB spectrum 
to constrain the 6 Li abundance produced via Galactic cosmic rays. By comparing this Galactic 
6 Li component to the observed solar abundance we can then place an upper limit on the residual 
6 Li which (presumably) was produced by SFCR. As described in §3.1, with the observed EGRB 
spectrum in hand we can place an upper limit on its fraction of pionic origin. In the case of 
SFCR-produced pionic gamma-rays, we can place constraints directly only in the presence of a 
model for the SFCR redshift history. Since a full model is unavailable, below (§5) we adopt the 
"single-redshift approximation." However, in the case of galactic cosmic rays we have a better 
understanding of the redshift history of the sources. Therefore, we will follow Pavlidou & Fields 
(2002) to calculate the pionic differential gamma-ray intensity for some set of energies 



/ 7 _ B = ■ ~ . / dz 



p*(z)L^[(l + z)E] 



AttHq^mw Jo y/n A + n M {i + z) 3 

( dt I dz) hA z) 

(26) 



1 ( 1 ^ £ dz(dt/dz)p*(z) 



Mo,mw VW),MW / J dz{dt/dz)p\{z) 

where I 7ir is in units of s _1 cm~ 2 GeV -1 , and Vmw is the present Milky Way star formation 
rate. For the pionic gamma-ray luminosity L 7ir we will, as before, use the pionic gamma-ray source 
function adopted from Pfrommer & Enfilin (2003) (a 7 = 2.75 for GCR spectrum), however we 
will let the normalization be determined by maximizing the pionic contribution to the EGRB. The 
adopted parameters, cosmology, and cosmic star formation rate we keep the same as in previous 
subsection. 

Once we obtain the spectrum we can then fit it with 

ln(I Jjr E 2 ) = -14.171 - 0.546 In E - 0.131(ln E) 2 + 0.032(ln Ef (27) 
where E is in GeV and / is in photons cm~ 2 s _1 sr _1 GcV -1 . The free leading term in the above 
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Fig. 1. — In the upper panel of this figure, we plot the pionic (dashed lines:green- maximized, blue- 
normalized to the Milky Way) EGRB spectrum, where decaying pions are of GCR origin, compared 
to the observed EGRB spectrum (solid line, fit to data); for purposes of illustration, we use the 
Strong, Moskalenko, & Reimer (2004) data points, which are given in red crosses. The bottom 
panel represents the residual function, that is, log[(7 E 2 ) b s / (J E 2 )^] = log(iobs/-^r)- 

equation is set by requiring that / 7ir = I 7 Q b S at the energy E = 0.44 GeV which maximizes pionic 
contribution by demanding that the pionic gamma-ray spectrum always stays below the observed 
one (since the feature of pionic peak is not observed). We also fit the Strong, Moskalenko, & Reimer 
(2004) data with 

ln(/ 7j0bs £ 2 ) = -14.003 - 0.144 In E - 0.097(ln£) 2 + 0.017(ln£) 3 (28) 
in the same units. 

By going through the procedure described in §3.1 (see Prodanovic & Fields 2004, for more 
detail) we can now obtain an upper limit to the fraction of pionic gamma-ray compared to the 
Strong, Moskalenko, & Reimer (2004) observed EGRB spectrum. This maximized pionic (green 
dashed line), as well as the observed, gamma-ray spectrum is presented in Fig. 1. We find the 
upper limit to pionic fraction to be f w (> O.lGeV) = f 01 dE J 77r / J Q 1 dE J 7)0 b S = 0.75. We note in 
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passing that a maximal pionic fraction as appears in Fig. 1 gives a poor fit at energies both above 
and below the matching near 0.4 GeV, suggesting the presence of other source mechanisms. This 
mismatch reflects a similar problem in the underlying Galactic 7-ray spectrum, and suggests that 
the pionic contribution to the EGRB is in fact sub-maximal. 



Thus, the pionic gamma-ray flux above 0.1 GeV is L /7T (> 0.1 GeV) = 0.83 
From eq. (21) it now follows that the total flux is I l7T (> 0) = 2.31 x 10~ 5 cm~ L s 
we can now use eq. (14) to find the GCR 6 Li mole fraction 



10 5 cm 2 s 1 sr 1 . 
sr^ 1 . As before, 



Y 6 



Li 



M>0) 



Yq t ■ I 9.06 x 10 -5 photons cm 2 s 1 sr 1 

> hl e / GCR 



0.25 



(29) 



and thus, SFCR-produced 6 Li can be at most (neglecting the 6 Li depletion) the residual 6 Li 

( y «Li \ 



V 6 



Lip 



o / SFCR 




Ll e / GCR 



(30) 



With the appropriate scaling between 7 Li and 6 Li as given in Table 1, we can then determine the 
total elemental Li = 7 Li + 6 Li abundance and compare it to the primordial values from (19) and 

(20): 

= 3.45 x 10~ 10 = 0.90 (-^) =2.81 (— M (31) 

SFCR V H / p,thy V H / p,obs 



H 



So far we have been determining the maximized pionic fraction of the EGRB based only 
on the shape of the pionic spectrum. However, in the case of normal galaxies we have a better 
understanding of what that fraction should be. That is, we can normalize pionic spectrum to 
the one of the Milky Way, and then integrate over the redshift history of sources. Following 
Pavlidou & Fields (2002) (and references therein) we set up the normalization by requiring that 
Jo 1 GeV dE W ( z = °) E ) = Jo 1 GeV rf %,Mw(£) = 2.85 x 10 42 s^ 1 . With the analytic fit of the 
shape of the pionic spectrum from Pfrommer & Enfilin (2003) this now gives: 



L~ (z = 0,E) = 9.52 x 10 44 s" 1 GeV^ 1 ( — \ ' + ( — \ 



- Oi'y 1 5*y 



(32) 



where a 7 = 2.75 for GCR spectrum, and <5 7 = 0.14a~ 1 ' 6 + 0.44. Now we can use equation (26) to 
obtain the pionic spectrum which is plotted on the Fig. l(blue dashed line) .We use star formation 
rate ipMW = 3.2M yr _1 (McKee 1989). Finally we find that in this case, when pionic spectrum is 
normalized to the Milky Way, the GCR 6 Li mole fraction that accompanies it is 



Ye 



Li 



Ye 



0.14 



Li 



(33) 



© / GCR 
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which then gives (Li/H) SF CR = 3.96 x 10~ 10 = 1.03( 7 Li/H) Pithy = 3.22( 7 Li/H) Pi0bs , which is of 
course a weaker limit than the maximal pionic case. 

We thus see that in a completely model-independent analysis, current observations allow the 
possibility that SFCRs are quite a significant source of 6 Li and of 7-rays. Indeed, we cannot 
exclude that SFCR-produced lithium can be a potentially large contaminant of the pre-Galactic Li 
component of halo stars, which would exacerbate the already troublesome disagreement with CMB- 
based estimates of primordial 7 Li. Consequently, we conclude that models for SFCR acceleration 
and propagation should include both 7-ray and 6 Li production; and more constraints on SFCR, 
both theoretically (e.g., space and time histories) and observationally (e.g., EGRB and possibly 
diffuse synchrotron measurement), will clarify the picture we have sketched. 

Note that had we also considered the possibility of depletion of 6 Li , we would have found a 
greater 6 Li residual, and thus had an even larger SFCR-produced component. 



In this section we turn to the as-yet unobserved cosmological component of cosmic rays, and 
to the synthesis of lithium by SFCR. This lithium component would be the first made after big 
bang nucleosynthesis. Any Li which is produced this way prior to the most metal-poor halo stars 
would amount to a pre-Galactic Li enrichment and thus would be a non-primordial Li component, 
unaccompanied by beryllium and boron production. This structure-formation Li would be an ad- 
ditional "contaminant" to the usual components in halo stars, the 7 Li abundance due to primordial 
nucleosynthesis, and the 6 Li and 7 Li contribution due to Galactic cosmic rays Ryan et al. (2000). 
Moreover, the pre-Galactic but non-primordial component would by itself be indistinguishable from 
the true primordial component, and thus would lead to an overestimate of the BBN 7 Li production. 

Our goal in this section is to exploit the 7-ray connection to constrain the structure-formation 
Li contamination. Unfortunately, we currently lack a detailed understanding of the amount and 
time-history of the structure formation cosmic rays (and resulting 7-rays and Li) . Thus we will make 
the conservative assumption that all structure formation cosmic rays, and the resulting 7s and Li, 
are generated prior to any halo stars. Furthermore, we will assume that the pionic contribution to 
the EGRB is entirely due to structure formation cosmic rays. This allows us to relate observational 
limits on the pionic EGRB to pre-Galactic Li. 

With this assumption and a SFCR composition « = 0.1, we can now use the 

appropriate scaling factor from Table 1 to rewrite eq. (13) 



5. Li and Gamma-Rays From Cosmological Cosmic Rays 




(34) 



(35) 
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(%)-^( ^tr-tV. ) 

where we used the solar lithium mole fraction Y( 6 Li) Q = 1.09 x 10~ 10 . 

To set up an extreme upper limit on pre-Galactic SFCR 6 Li, we assume that the entire pionic 
extragalactic gamma-ray background came from SFCR-made pions, and was created prior to any 
halo star. As mentioned in the previous section, the method used in subtraction of the Galactic 
foreground is crucial for obtaining the EGRB spectrum. What is more, the EGRB spectrum is 
an important input parameter in the Prodanovic k, Fields (2004) analysis whose estimates of the 
maximal pionic gamma-ray flux we will use here. Our results for the SFCR lithium upper limits 
are collected in Table 2. The results depend on the choice of the EGRB spectrum as well as the 
redshift of origin of cosmic-rays according to the single-redshift approximation used by Prodanovic 
& Fields (2004) to obtain the maximal pionic EGRB fraction. Note that we considered only the 
two most extreme redshifts to illustrate the results. In the Table 2, z is the redshift, i 7 , 7r (> 0) is 
the upper limit for the pionic 7-ray intensity above energy determined from (18) as explained in 
the previous section, (Li/H)g^ R is the upper limit to total ( 6 Li+ Li) lithium abundance that can 
be of SFCR origin, while Li* hco and Li° bs are the theoretical and observational primordial lithium 
abundances respectively as given in equations (20) and (19). 

Notice that for the case of Keshet, Waxman, & Loeb (2003) EGRB, since a spectrum was 
unavailable, the procedure described in the section §3.1 for maximizing the pionic fraction of the 
EGRB could not be used. Thus, to place an upper limit on SFCR lithium we assumed that the 
entire EGRB can be attributed to decays of pions, that is, assume l 1 = I 7)Vr o. For the Sreekumar 
et al. (1998) and Strong, Moskalenko, & Reimer (2004) EGRB spectra, we use the upper limits to 
I 7j7r o obtained by Prodanovic & Fields (2004). Once the I 7)7r o is set we can use (36) to find the 
SFCR 6 Li upper limit. 

To find the total halo star contribution we must also include 7 Li, which is in fact produced 
more than 6 Li in aa fusion: as seen in Table 1, ( 7 Li/ 6 Li) gFCR = (v??-) / (0%? .) ~ 2. The total 
SFCR elemental Li production appears in Table 2, both in terms of the absolute Li/H abundance 
and its ratio to the different measures of primordial Li (§3.2). 

From Table 2 we see that the maximal possible SFCR contribution to halo star lithium could be 
quite substantial. If the pre-Galactic SFCR component is dominantly produced at high redshift (i.e., 
as in the z ~ 10 results) then the maximum allowed Li production is can exceed the primordial Li 
production (however it is estimated), in some cases by a factor up to 25! The situation is somewhat 
better if the pre-Galactic SFCR production is at low redshift, but here it is hard to understand 
how this would predate the halo star component of our Galaxy. The high-redshift result is thus 
the more likely one, but also somewhat troubling in that the limit is not constraining. The indirect 
limits on SFCR Li in the previous section are somewhat stronger, but these also hold the door open 
for a significant level of pre-Galactic synthesis. 
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Table 2: Upper limit on Li of SFCR origin 



EGRB [cm- 2 s- 1 sr- 1 ] 


z 


hA> °) 


(Li/H)sFCR 


— 1 

SFCR 
T jthco 


— 1 ;iiiax 

SFCR 
Lig bs 


u 


Sreekumar et al. (1998) 
/ 7 , obs (> 0.1) = 1.57 x 10- 5 



10 


8.78 x 10~ 6 
1.22 x 10" 4 


2.19 x 10" 10 
3.04 x 10~ 9 


0.57 
7.95 


1.78 
24.7 


0.91 
0.15 


Strong et al. (2004) 
/ 7 , obs (> 0.1) = 1.11 x KT 5 




10 


4.59 x 10~ 6 
6.27 x 10~ 5 


1.14 x 10~ 10 
1.56 x 10~ 9 


0.30 
4.09 


0.93 
12.69 


1.29 
0.21 


Keshet et al. (2003) 
/ 7 ,obs(> 0.1) < 0.5 x 10~ 5 




10 


< 6.5 x HT 6 

< 4.03 x 10~ 5 


< 1.62 x 10" 10 

< 1.00 x 10~ 9 


0.42 
2.63 


1.32 
8.16 


2.86 
0.46 



We caution that the lack of a strong constraint on SFCR Li production is not the same as 
positive evidence that the production was large. Recall that we have made several assumptions 
which purposely maximize the SFCR contribution; to the extent that these assumptions fail, the 
contribution falls, perhaps drastically. A more detailed theoretical and observational understanding 
of the SFCR history, and of the EGRB, will help to clarify this situation. Moreover, given that the 
halo star Li is already found to be below the CMB-based 7 Li BBN results, we are already strongly 
biased to believe that the pre-Galactic SFCR component is not very large. Thus one might be 
tempted instead to go the other way and use Li abundances to constrain SFCR activity. 

We thus now go the other way and use solar 6 Li to constrain the SFCR 7-ray flux. Again, given 
our incomplete knowledge of SFCRs, we must adopt a simplifying assumption about the degree 
of 6 Li production which is due to SFCR. To be conservative, we make the extreme assumption is 
that all of the solar 6 Li is produced by SFCR, and thus find via eq. (34) that 7-flux is J 7 7r o(> 
GeV) > 1.86 x 10~ 5 photons cm~ 2 s _1 sr _1 . From (18) we can determine I 7j7r o(> 0.1 GeV) to 
be 0.23 x 10~ 5 < i" 7)7r o(> 0.1 GeV) < 1.43 x 10" 5 photons cm~ 2 s _1 sr _1 depending on the redshift 
of pionic 7-rays, which is below the observed level as determined by Sreekumar et al. (1998), and 
a factor of 2-14 lower than the prediction based on GCR. Thus, for a given observed intensity 
I° bs (> eo) we can now u se (18) to constraint the hadronic fraction of EGRB, that is, calculate 
fir(> e o, z ) which is also presented in the Table 2. 

However, since Li is being depleted, the use of the solar 6 Li abundance does not give us 
the upper most limit to the required pionic gamma-ray flux 7 7 7r o(> 0). Thus, if one would to 
compensate for the depletion, the pionic fraction f n (> eo, z) would become even larger. 

Indeed, this may suggest a solution to the EGRB overproduction by GCRs, seen in the previous 
section. If 6 Li is mostly made by SFCRs, then the associated 7-ray production is in line with the 
observed background. In this case, 6 Li would still be of cosmic-ray origin, but not dominated by 
GCR production. Such a scenario faces tests regarding 6 Li and other LiBeB abundances and their 
Galactic evolution. A detailed discussion of this scenario will appear in a forthcoming study. 
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6. 



Discussion 



The main result of this paper is to identify and quantify the tight connection between "Li and 
the EGRB as measures of cosmic-ray history. Specifically, these two observables provide measures 
of average, gas-weighted cosmic-ray fluence. Moreover, the observables are complementary, in that 
6 Li samples local fluence, while the EGRB encodes the cosmic mean fluence. 

We present scaling laws which relate 6 Li and the EGRB intensity for different cosmic ray 
spectra appropriate for GCR and SFCR populations. Using these scalings, and assuming that our 
local 6 Li measurements are typical, we can test the self-consistency of 6 Li and EGRB observations 
in a relatively model-independent manner. We find that if SFCR dominate the pionic EGRB, then 
the associated 6 Li production can be a significant and perhaps dominant contribution to the solar 
abundance. On the other hand, we find that if 6 Li production is dominated by GCRs, then the 
associated 7-ray production is enormous, at least a factor of two above the observed intensity. 

Furthermore, using the EGRB we use two different lines of argument to place an upper limit 
on the SFCR contribution to pre-Galactic lithium in halo stars. Such a component of lithium would 
be confused with the true primordial abundance and thus would exacerbate the existing deficit in 
halo star Li relative to the CMB-based expectations of BBN theory. Unfortunately, current EGRB 
data are such that our model-independent upper limit (which most assume, among other things, 
that all SFCRs are created prior to any halo stars) is very weak. In particular, we cannot exclude 
the possibility that a significant portion of pre-Galactic lithium is due to SFCRs. We thus find that 
the nucleosynthesis aspects of SFCRs are important and deserver further more detailed study. 

A full understanding of the implications of the relationships among 6 Li, diffuse 7-rays, and 
cosmic ray populations thus awaits better observational constraints (both light elements and espe- 
cially the EGRB) as well as a more detailed study of SFCRs. Having shown the importance of both 
the 6 Li and 7-ray constraints, it is our hope that these observables will be calculated in models of 
galactic and structure formation cosmic rays, and that both 6 Li and the EGRB are used in concert 
to constrain cosmic-ray interactions with diffuse matter. The results of these models will go far to 
address some of the questions which this study has raised. 

We are grateful for illuminating discussions with Vasiliki Pavlidou. This material is based 
upon work supported by the National Science Foundation under Grant No. AST-0092939. 



The interactions of cosmic-ray species i with target nucleus j produces species A; at a rate per 
target particle of 



A. 



Notation and Normalization Conventions 




(Al) 
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Here E is the cosmic-ray energy per nucleoli, (Jij^k is the energy-dependent production cross 
section, with threshold £th,fc) and fa is the cosmic-ray flux. The rate per unit volume for i + j — > k 
is thus qk = TkUj. 

Note that the flux in eq. (Al) is position- and time-dependent. To isolate this dependence, it 
is useful to define a total, energy-integrated, flux 

$i= I dEfa (A2) 

where we choose the lower integration limit to always be the minimum threshold £th,min for all 
reactions considered; in our case this is the a + a — ► 7 Li threshold of 8.7 MeV/nucleon. From 
eqs. (Al) and (A2) it follows that 

Oij^k = Tk/^i (A3) 

represents a flux-averaged cross section. Also note that if the spectral shape of fa is constant (as 
we always assume), then so is dij-^k, an d the flux $j contains all of the time and space variation 

of rv 

Finally, two conventions are useful for quantifying abundances. Species i, with number density 
rii, has a "mole fraction" (or baryon fraction) Y{ = ni/n^. It is also convenient to introduce the 
"hydrogen ratio" = nj/n H = Yi/Y R . 



B. Cosmic Gamma-Radiation Transfer 

The expression for 7-ray intensity in a Friedmann universe is well-known (Stecker 1969), but 
usually expressed in redshift space. For our purposes, the result expressed in the time domain is 
critical, and indeed is more fundamental, so we give a derivation based on the Boltzmann equation. 
For this section we adopt units in which c = 1. 

The differential photon (number) intensity I is directly related via 

I(p,x,t)=p 2 f(p,x,t) (Bl) 

to the 7-ray distribution function f(p,x,t) = d?N/d 3 pd 3 x Here p and x, as well as the volume ele- 
ments, are physical quantities (and thus subject to change with cosmic expansion). The distribution 
function is related to the photon sources via the relativistic Boltzmann equation 

p»d,f-Tf P yPf = E^L (B2) 

where gravitational effects enter through the Affine connection T, where E = p = \p\, and where 
the source function (number of photons created per unit volume per unit time) is q. 

For an isotropic FRW universe we have / = f(E,t), and thus 
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where q(E) = dq/dE and where we neglect photon scattering and absorption. 

We now note that a given photon's energy E drops due to redshifting as a -1 . It is thus useful 
to define a comoving energy e = aE; with a(to) = 1, we see that e is also the present-day (observed) 
photon energy. Changing variables from f(E, t) to /(e, t), and similarly for q, the energy-dependent 
d e term drops out; this is physically reasonable since we do not allow for scattering processes, and 
thus a photon's energy can only change due to redshifting. We then have 

9tf - a "4^" (B4) 
which, for any fixed comoving energy e, integrates to 

f(e,t) = ^J*dt'a 2 q(e/a) (B5) 
Equation (Bl) then gives the intensity 

where q com = a 3 q is the comoving source rate. Equation (B6) is the usual expression (which is 
often then expressed in terms of an integral over redshift). Finally, if we integrate over the entire 
energy spectrum, and evaluate at the present epoch to (when ao = 1), we have 

/(>0,t )= / del(e,t ) = — dt'q com (>0) (B7) 
Jo 47r Jo 

where q CO m(> 0) = f °°du q{u) is the total source rate, integrated over rest-frame energy. 

We see from eq. (B7) that the energy-integrated intensity is the same as one would find from 
uniform sources in a non-expanding universe (which have been "switched on" for a duration t). 
This result is physically sensible, because the two effects of cosmic expansion are to introduce a 
particle horizon and redshifting. The energy integration removes the effect of redshifting, so that 
the only effect is that of the particle horizon, which acts to set the integration timescale. 
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